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Nonlinear Spinor Fields in Bianchi type-VI space-time 

Bijan Saha 

Laboratory of Information Technologies 
Joint Institute for Nuclear Research 
141980 Dubna, Moscow region, 

Within the scope of Bianchi type-VI cosmological model we study the role of spinor 
field in the evolution of the Universe. It is found that due to the spinor affine connections 
the energy-momentum tensor of the spinor field possesses non-diagonal components. The 
non-triviality of non-diagonal components of the energy-momentum tensor imposes some 
severe restrictions either on the spinor field or on the metric functions or on both of them. 
But unlike in cases of Bianchi type-I or VIq, in case of Bianchi type-VI model it does not lead 
to the elimination of spinor field nonlinearity or mass term in the spinor field Lagrangian. It 
is also found that depending on the sign of self-coupling constant the model can give rise to 
late time acceleration or generate oscillatory mode of evolution. 

PACS numbers: 98.80.Cq 

Keywords: Spinor field, late time acceleration, oscillatory solution, anisotropic cosmological mod¬ 
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I. INTRODUCTION 


With the more and more observational data available from far sky, the need for a change in the 
standard cosmological paradigm becomes inevitable. Prior to 1998 when we had no idea about 
the accelerating mode of expansion the available observational data were well fit in the standard 
Einstein model. But the discovery and further reconfirmation of the existence of the late time 
accelerated mode of expansion [1,2] have opened a new window for change. Along with that 
comes out a number of altemativ; models of the evolution of the Universe. 

The most popular among the models are those which consider the old Einstein theory with a 
new ’’matter” as a source field. The models with A-term [3-51, quintessence [6-111. Chaplygin 
gas r 12-191 etc. are among the most studied ones, though sojrr e other models of daik energy are 
also proposed. After some remarkable works by different au ijicrs [20—34], showing the important 
role of s pir or field in the evolution of the Universe, it has been extensively used to model the dark 
energy. This success is directly related to its ability to answer some fundamental questions of mod¬ 
ern cosmology: (i) Problem of initial singularity and its possible elimination 122-26, 35-381: (ii) 


problem of isotropization [24, 25, 27, 36 
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'rom ekpyrotic matter to phantom matter, as well as Chaplygin 
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re noticed that in earlier works only the diagonal components of the energy- 
tensor of the spinor field were taken into account. But recently it was shown that 
due to its specific behavior in curve spacetime the spinor field can significantly change not only 
the geometry of spacetime but itself as well. The existence of nontrivial non-diagonal components 
of the energy-momentum tensor plays a vital role in this matter. In [47, 48] it was shown that de¬ 


pending on the type restriction imposed on the non-diagonal components of the energy-momentum 
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tensor, the initially Bianchi type-I evolves into a LRS Bianchi type-I spacetime or FRW one from 
the very beginning, whereas the model may describe a general Bianchi type-I spacetime but in that 
case the spinor field becomes massless and linear. The same thing happens for a Bianchi type- 
Vto spacetime, i.e., the geometry of Bianchi type-V/o spacetime does not allow the existence of a 
massive and/or nonlinear spinor field [49]. 

Anisotropic Bianchi type VI cosmological models were studied by many authors [50-54]. In 


this report we study the role of spinor 
cosmological model. 


field in the evolution of a Bianchi type VI an 


sotropic 


II. BASIC EQUATION 


Let us consider the case when the anisotropic space-time is filled with nonlinear spinor field. 
The corresponding action can be given by 


y{g\ y, \jr) = J L^idCl 

with 

L = Lg T L sp . 

Here L g corresponds to the gravitational field 

L -A 

g- 2 k’ 


( 2 . 1 ) 

( 2 . 2 ) 


(2.3) 


where R is the scalar curvature, K — 8 kG, with G being Einstein’s gravitational constant and L sp 
is the spinor field Lagrangian. 


A. Gravitational field 


The gravitational field in our case is given by a Bianchi type-VI anisotropic space time: 

ds 2 = dr — a\e 2mX3 dx 2 — a\e lnx ' i dx\ — a\dx\, (2.4) 

with ai, «2 and <23 being the functions of t ime only and m and n are some arbitrary constants. 

The nontrivial Christoffel symbols for (12.41) are 
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-m, r 32 = n, 
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= a2d2e lnx 3 , r^ 3 = a 3 a 3 , 

3 _ma 2 2mX3 r 3 
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(2.5) 
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The nonzero components of the Einstein tensor corresponding to the metric (12.41) are 

_ _<h _03 _ 0203 

1 02 O3 02 03 a 2 ’ 

2 03 o\ a,3 a\ m 2 

&2 — I 2> 

03 a i a 3 a \ a 3 

3 a | 02 a i ci2 mn 

C?3 2 “ ■ 

a i «2 ai «2 <33 

g dj ^2 «2 ^3 ^3 di m 2 —mn+n 2 

G 0 — -1- 2 -> 

a\ 02 02 03 03 a\ a 3 


G 


0 

3 


(m 


d3 di d2 

n) - m - 1 - n —. 

03 a\ 02 


(2.6a) 

(2.6b) 

(2.6c) 

(2.6d) 

(2.6e) 


B. Spinor field 


For a spinor field t ff, the symmetry between y/ and y/ appears to demand that one should choose 
the symmetrized Lagrangian (55]. Keeping this in mind we choose the spinor field Lagrangian as 
[24]: 


-sp 


LA y\ 


m 


sp 


W-F, 


(2.7) 


where the nonlinear term F describes the self-interaction of a spinor field and can be presented 
as some arbitrary functions of invariants generated from the real bilinear forms of a spinor field. 
Since yr and y/* (complex conjugate of yr) have four component each, one can construct 4x4=16 
independent bilinear combinations. They are 


S — iffy/ (scalar), (2.8a) 

P = lyrj’yf (pseudoscalar), (2.8b) 

= (y/y^y/j (vector), (2.8c) 

= (y/y^yr) (pseudovector), (2.8d) 

Q^ v — (y/o^ v yf) (antisymmetric tensor), (2.8e) 


where a^ v — (z/2)[y^y v — y v y u j. Invariants, corresponding to the bilinear forms, are 


/ = S 2 , (2.9a) 

J = P 2 , (2.9b) 

I v = v M v^ = (y/y^ y/j gji v ( YY V Y) , ( 2 . 9 c) 

I A = A m A^ = (vr’fw)gnv(wr^7 v v); (2.9d) 

Iq = Q^vQ^ = (y/r7 MV y/)^ a g Vj6 (yAf7“ /3 yA). (2.9e) 


According to the Fierz identity, among the five invariants only I and J are independent as 
all others can be expressed by them: I v — —I A — I + J and Iq = I — J. Therefore, we choose 
the nonlinear term F to be the function of I and J only, i.e., F = F (/, J), thus claiming that it 
describes the nonlinearity in its most general form. Indeed, without losing generality we can 
choose F = F(K), with K taking one of the following expressions {/, J, I + J, I — i}.. Here is 
the covariant derivative of spinor field: 
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with being the spinor affine connection. In (12.71) y’s are the Dirac matrices in curve space-time 
and obey the following algebra 

y^y v + 7 V 7 /J = 2g MV (2.11) 

and are connected with the flat space-time Dirac matrices y in the following way 

8fiv(x) = e a fJL {x)e b v (x)rj ab - 7m W = (2.12) 

where e" is a set of tetrad 4-vectors. 

For the metric (12.41) we choose the tetrad as follows: 

= 1, e^—a\e~ mX3 : = a 3 . (2.13) 

The Dirac matrices y^(x) of Bianchi type-VI space-time are connected with those of 
Minkowski one as follows: 


with 


jnx 3 —nx 3 i 

7° = 7°, 7 1 = -?\ T 2 =-T 2 ) T 3 = —T 3 

a\ ci2 a 3 

T 5 = -«y = i’7 D 7 1 7 2 7 3 = -iffff = T 5 


7° 



f 


0 o' \ 

—o' 0 J’ 


where o, are the Pauli matrices: 





Note that the y and the o matrices obey the following properties: 

??' + r'Y = 2r] ij ; i,j = 0 ,1,2,3 
Yf + fY = 0, (y 5 ) 2 = /, i = 0, 1,2,3 

o J o k = 8 jk + i£ jk ,o l , j 1 k^l = 1,2,3 

where i],/ = {1,-1, —1,-1} is the diagonal matrix, 8j k is the Kronekar symbol and £ jk i is the 
totally antisymmetric matrix with £123 = + 1 . 

The spinor affine connection matrices F^, (x) are uniquely determined up to an additive multiple 
of the unit matrix by the equation 


0 -r^yp -r M y v + y v r M = 0, (2.14) 

with the solution 

r„ = Ijir'Vv 1 - \i P f rjv (2.15) 

From the Bianchi type-VI metric (12.151) one finds the following expressions for spinor affine 
connections: 


r 0 = 0, (2.16a) 

ri = ~(d\— e~ mX3 , (2.16b) 

2 V a 3 J 

r 2 - + (2.16c) 

r 3 = y 7 3 f- (2.16d) 
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C. Field equations 

Variation of (12.11) with respect to the metric function g^ v gives the Einstein field equation 

g;=k;-^s;r=-kt*, an > 

where R ^ and R are the Ricci tensor and Ricci scalar, respectively. Here 7’ u v is the energy momen¬ 
tum tensor of th e spinor field. 

Varying (12.71) with respect to xfr( y/) one finds the spinor field equations: 

?y /i V jU t//-m sp v/-^t/A-j^y 5 i/A = 0, (2.18a) 

i^gYY P + m s,pY + @y + &yy 5 = 0, (2.18b) 

where we denote & — 2SFkKj and { F — 2PFkKj, with Fk — dF/dK, Kj = dK / dl and Kj = dK / d.J. 
In view of (12.181) can be rewritten as 

U P = ^[V7fVpV-VpVy ll v]~m sp W-F(I,J) 

= ^Vlr^nV-m sp Y] - \\7nvf + m sp y]¥-F(I,J), 

= 2{IF i + JFj)-F = 2KF k -F{K), (2.19) 

where K — {/, 7,7 + 7,7-7}. 


D. Energy momentum tensor of the spinor field 

The energy-momentum tensor of the spinor field is given by 

T p = ^g pv {YYgVvY+WvVgY-VgYYvY-VvYYgY) ~^L sp . (2.20) 

Then in view of (12.101) and (12.191) the energy-momentum tensor of the spinor field can be written 
as 

T/ = ^g pv (YYgdvY+YYvdgY-d^YYvY-dvYYgY) 

- ^g pv v{Y^v +r v r M + Yv^n +r M rv) w-s p ( 2 kf k - f(k)) . ( 2 . 21 ) 

As is seen from (12.211) . is case if for a given metric F u \s are different, there arise nontrivial non¬ 
diagonal components of the energy momentum tensor. 

We consider the case when the spinor field depends on t only, i.e. y/ = y/(/). Then inserting 
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(12.101) into (12.211) one finds 

To = ^g°° {Y7 oY~ ^YoY)-L sp, (2.22a) 

T\ = -^“^(rin+rirOvz-Lsp, ( 2 . 22 b) 

T 2 = -^g 22 Y(y2T2 + r 2 Y2)Y-U pi (2.22c) 

T 3 3 = ~g 33 ¥( 73^3 + ^ 373 ) ¥-L sp , (2.22d) 

T? = \g °° {WY\W- fY\ Y) ~ ^g°°W(YoT\ + H 70) y, ( 2 . 226 ) 

t 2 = \g 00 (vy 2 ¥-^72 y) - ^g 00 y (ror 2 +r 2 ro) y, (2.22f) 

r 3 ° = ^ 00 (y73Y-Y73y) - ^g 00 Y (ror 3 +r 3 ro) y, (2.22 g ) 

T\ — ~-^g li Y(Y2Ti +Ti72 + 7 ^ 2 + T27i) Y> (2.22h) 

r 3 2 = -^ 22 t/A(73r2+r2r 3 + 72r 3 +r 3 7 2 )t/A, (2.22i) 

73 = -^g 1 V(73ri+r 1 7 3 + 7 1 r 3 +r 3 7i)t/A. (2.22j) 


Further inserting (12.161) into (12.221) after a little manipulations for the components of the energy- 
momentum tensor one finds: 

yO _ 

— 

7 )' = 

if = 

j^O _ 

2 — 

T’O _ 

'3 — 

T 1 — 
i 2 — 


T 1 — 

1 3 — 

'3 — 

As one sees from (12.221) and (12.231) the non-triviality of non-diagonal components of the energy 
momentum tensors is directly connected with the affine spinor connections I/S. 


m sp S + F(K), 

T 2 = T 3 = F(K) — 2 KF k , 

Yf 3 f 1 y°Y — 


me "“ 3 ci 1 


ne mX3 a\ 2 


«3 


-A z , 


ime nx3 «2 


4 a 3 

0 , 

|e (m+n)*a a2 

4 ai 

e ( m +«)-T3 a2 

4 ai 
a 3 / a 3 


Y7 2 7 3p / ) Y — 


m e' m a 2 j 


4 a 3 


03 

A 1 , 


---) YfffY- 

a 1 a 2 ) 

£1 _ £2 \ A 3 _ m + w A o~ 


m + n 


CI 3 


Y7 l 7 2 7 3 Y 


a 1 «2 

di 


«3 


4 ai \a 3 ai 


e " U ' 3 CI 3 


1 e nX3 CI 3 / a 2 


a 2 \a 2 


£3 

03 


V^T 2 T 3 7° = 


£3 _ £1 \ a 2 
4 aj ^a 3 ai/ 




e nX3 a 3 / a 2 a 3 


4 a 2 \ci 2 C13 J 


A 1 . 


(2.23a) 

(2.23b) 

(2.23c) 

(2.23d) 

(2.23e) 


(2.23f) 

(2.23g) 

(2.23h) 
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From (12.181) one can write the equations for bilinear spinor forms (12.81) : 


,4° 


■3 


m — n 


a 3 


we also introduce the volume scale 


S 0 + &A ° 0 = 0, 

(2.24a) 

fl)-<J>A§ = 0, 

(2.24b) 

f — n -> 

-- A 3 0 + <S>P 0 -&S 0 = 0, 

«3 

(2.24c) 

. 0 m — n , n 

A5 A 0 — 0, 

<33 

(2.24d) 

n m — n 3 

Vo- v 3 0 = 0, 

<33 

(2.24e) 

-v°o + ^Ql 0 +^Ql l = 0, 

(2.24f) 

e 30 -d>v 3 = 0, 

(2.24g) 

= 0, 

(2.24h) 

= A^V, v% = v»V, Q% v = Q^ v V and <S> = 

m s p + @. Here 

V = a\a 2 ci 3 . 

(2.25) 

and taking the first integral one gets 


+ (A°) 2 -(A 3 ) 2 = Ci = Const, 

(2.26a) 

2 + ( v o) 2 -( v o) 2 = C 2 = Const 

(2.26b) 


(: Qo) 2 + (Ql l 

Now let us consider the Einstein field equations. In view of (12.61) and (12.231) we find the follow¬ 
ing system of Einstein Equations 


a 2 ^ ®3 <3 2 a 3 
Cl2 CI3 Cl 2 0-3 
CI3 fli CI3 Cl l 

<23 a\ CI 3 ci\ 


m 


k(F(K)- 2KF k ), 


= k{F{K)- 2KF k ), 


a\ a 2 a, a 2 mn 

-1-1- h—y = k(F(K)- 2KF k ), 

a 1 a 2 a 1 a 2 ai 


a 1 a 2 a 2 0,3 CI 3 a 1 
ci\ ci 2 ci 2 (13 0,3 a 1 


hi — mn + n~ 




= K (m sp S + F(K)) , 


, x a 3 a\ a 2 

(m — n )- m - Yn — = 0 , 

CI3 a\ a 2 


0 = 
0 = 

0 = 
0 = 
0 = 


ne~ mx i a 1 


m e 


nx 3 


«3 

—A 1 , 


4 <23 

e (m+n)x 3 a2 


e mx ^ CI 3 


4 

e ~ nx 3 


a 1 
03 

a\ \a 2 


CI 3 


4 a 2 


a 2 

a 2 


( d\ 

<32 \ 


<32/ 

«i\ 

A 2 , 

ai) 


CI3 

V 


o m + n a 

A 3 -'—A 0 

a 3 


CI 3 


(2.27a) 
(2.27b) 
(2.27c) 
(2.27d) 
(2.21 o) 
(2.27f) 
(2.27g) 
(2.27h) 
(2.27i) 
(2.27j) 
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From (I2.27fl) and (|2.27g[) one dully finds 


A 2 = 0, and A 1 = 0. 


(2.28) 


In view of (12.281) the relations (I2.27il) and (|2.27j[) fulfill even without imposing restrictions on the 
metric functions. From (I2.27hl) one finds the following relations between A 0 and A 3 : 


a i a 2 


\ A :3 = '» + » 


Inserting (12.291) into (I2.24dl) one finds 


m + n A^ 


777 


-77 A 3 


771 
Cl 1 


03 


(72 

(7 2 


(2.29) 


(2.30) 


with the solution 


(Aq) = X 03 j , X 03 = const. (2.31) 

On the other hand from (I2.27el) one finds the following relation between the metric functions 


03 = -Xo 


1 / (m—r, 


(2.32) 


Thus the non-diagonal components of Einstein equations not only connected the different met¬ 
ric functions as was found in [261 . but also imposes some restrictions on the components of the 
spinor field. 

explicitly we have to address the diagonal components of Einstein 


To find the metric functions 


system. Explicit presence of <73 force us to impose some additional conditions. In an early work 
[261 we propose two different situations, namely, set (73 = y/V and 7/3 = V which allows us to 
obtain exact solutions for the metric functions. 

In a recent paper we imposed the proportionality condition, widely used in literature. Demand¬ 
ing that the expansion is proportion to a component of the shear tensor, namely 


&=N 3 o; 


(2.33) 


The motivation behind assuming this condition is explained with reference to Thome [561 . The 
observations of the velocity-red-shift relation for extragalactic sources suggest that Hubble ex¬ 
pansion of the universe is isotropic today within « 30 per cent [57. 58]. To put more precisely, 
red-shift studies place the limit 

(7 

— <0.3, 


H 


(2.34) 


on the ratio of shear <7 to Hubble constant H in the neighborhood of our Galaxy today. Collins 
et al. [59] have pointed out that for spatially homogeneous metric, the normal congruence to the 
homogeneous expansion satisfies the condition ^ is constant. Under this proportionality condition 
it was also found that the energy-momentum distribution of the model is strictly isotropic, which 
is absolutely true for our case. 

Let us now find expansion and shear for BVI metric. The expansion is given by 


= Uu + P 


lia 1 ' 


(2.35) 






and the shear is given by 


(2.36) 
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with 


°" v = 2 L 


Uu.aPv J r u v:aPu 


y'&PlJ.V, 


where the projection vector P: 


P P) P [IV SflV Py 8 y U^Uy. 

In comoving system we have = (1,0,0,0). In this case one finds 

„ a\ do do V 

& —-1-1-= 77; 

a\ a 2 a 3 V 

and 


(2.37) 


(2.38) 


(2.39) 


CTi = -- 


2 _ 1 
° 2 3 

ol = -l 


— 2 — 

+ 

1 

1 NJ 

+-) 

di 

-k 

(2.40a) 

a\ 

«2 

<33/ 


3 

—2— 

do 

+ — 

+ —) 

d 2 

-k 

(2.40b) 

Cl2 

ao, 

a\J 

a 2 

3 


_ do, 
— 2 —- 

d\ 

+ — 

+ 

1 

1 b° 

_ d 3 

-k 

(2.40c) 

«3 

ai 

a 2 / 

ao, 

3 


One then finds 


C7“ = 


3 

E 


i=\ \ a i 


-tT 


i= 1 




Inserting (12.321) into (12.391) . (12.401) and (12.411) we find 


2 m — nd\ m — 2n do 
& = -- +-- 


and 


= 

°2 = 


CJ 3 3 = 


m + n 


m 

— n ao 

a\ 

d 2 \ 

a\ 

a 2 J 

a\ 

d 2 \ 

a\ 

a 2 J 

a\ 

d 2 \ 

a\ 

a 2 J 


On account of (12.321) . (I2.40cl) . (12.251) from (12.331) one finds 


a i 


V 


—In 


_i_ |_at m—zn 

n—n )/ (m—2n) "1 3 /V3 m+n 

-v 


a 2 — X\ 


X 


X 0 


r _1_|_ tv j n—2m 

{?n—n)/(m—2n) "| 3 + “»+i 

-V 


«3 = 


*0 


—n/(m—n) 


X 


(m—ri)/{m—2n) 

1 


Xq 


■V 


Un 3 


(2.41) 

(2.42) 

(2.43a) 

(2.43b) 

(2.43c) 


(2.44a) 

(2.44b) 

(2.44c) 
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where X\ is an integration constant. Further taking into account that V = a i< 22^3 from (12.441) one 
finds 

m—n I m—2n 

X m-2n m—n = j (2.45) 

with either Xi = 1 or ^ 5 ^ + = 0. Since (m — n ) 2 + (m — 2 «) 2 ^ 0, we conclude Xi = 1. 

Hence for the metric functions finally we obtain 


(31 = 


V 

*o 


Un 3 


m—2n 
m+n 


Cl2 


V 

Xo 


k+N 3 


n—2m 
m+n 


a-i = Xq 


V 

*o 


Un 3 


(2.46) 


The equation for V can be found from the Einstein Equation (12.61) which for some manipulation 
looks 


V = 2 (m 2 - mn + n 2 )X™ 3 4 /V/ 3-2 *' 3 + [m sp S + 2 ( F(K ) - KF k )\ V. (2.47) 

In order to solve (12.471) we have to know the relation between the spinor and the gravitational 
fields. Let us first find those relations for different K. Let us recall that K takes one of the following 
expressions {/, /, / + 7, 1 — J}, with — 2 SF k K[ and & = 2 PF k Kj. 

In case of K — I, i.e. = 0 from (I2.24al) we duly have 

So = 0, (2.48) 

with the solution 

7 Vn y 0 

K = I = S 2 = =s> S = y, Vo = const. (2.49) 

In this case spinor field can be either massive or massless. 

As far as case with K taking one of the expressions {J, I+ J, I — 7} that gives Kj — ±1 is 
concerned, it can be solved exactly only fo r a ma ssless spinor field. 

In case of K — J, i.e. <f> = = 0 from (I2.24bl) we duly have 


with the solution 


Pq = 0 , 


7 Vn 2 V 0 

K = J = P 2 = => P = y 0 = const. 


In case of K = I + J the equations (I2.24al) and (I2.24bl) can be rewritten as 

S 0 + 2PF k A° 0 = 0, 


which can be rearranged as 


P 0 -2SF k A° 0 = 0, 


S 0 S 0 + PoPo = j f (55 + Eq 2 ) = j f ( V 2 K) = 0, 


with the solution 


V 2 

K = ^ 2 - Vo = const. 


Note that one can represent S and P as follows: 


c Vo Vo „ 

S = —sintk P — —cost). 

y ’ y 


(2.50) 

(2.51) 

(2.52a) 

(2.52b) 

(2.53) 

(2.54) 


(2.55) 
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The term 0 can be determined from (I2.52al) or (12.52bl) on account of (12.291) . (12.311) and (12.461) . 
Finally, for K — I — J the equations (I2.24al) and (12.2451) can be rewritten as 


S 0 -2PF K A° 0 = 0, 

Po-2SF k A 0 0 = 0, 

which can be rearranged as 

S 0 S 0 - P a P„ = 2 (sg _ />2) = £ (v 2 k) = 0, 


(2.56a) 

(2.56b) 


(2.57) 


with the solution 


K 


V 2 

v o 

V 2 ’ 


Vq — const. 


As in previous case one can rewrite S and P as follows: 


(2.58) 


S = —cosh0, P — —sinh0. (2.59) 

V V V 

Like previous case 0 can be determined from (I2.56al) or (I2.56bl) on account of (12.291) . (12.311) and 

C39 . 


III. SOLUTION TO THE FIELD EQUATIONS 


In thi s sec tion we sol ve the field equations. L et us b egin with the spinor field equations. In 
view of (12.101) and (12.161) the spinor field equation (12. 1 8al) takes the form 

if = 0, (3.1a) 

f ^ 1 V _\ n _ m — n 

As we have already mentioned, i/a is a function of t only. We consider the 4-component spinor 
field given by 


YJ 5 + ^Y + ^YY 5 = 0 . 


(3.1b) 


( Y l \ 
Yi 
V3 

W/ 


(3.2) 


Denoting </>,• = y/V Yi and Aq 


(m 


n)X^ 


from (13.11) for the spinor field we find we find 


01 -H<F0i + 

02 + l 9? 02 - 

03 — l T>03 + 
04 — l T>04 — 


2yi/3+A 3 


*0 


2V 1 /3-LV 3 


*0 


2U 1 /3+a 3 

*0 



1 2V'/ 3+/v 3 


03 

04 

01 

02 


0, 

0, 

0, 

0. 


(3.3a) 

(3.3b) 

(3.3c) 

(3.3d) 
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Further denoting 3/ = we can write the foregoing system of equation in the form: 


with 0 — col ( 0 1? 02 , 03 , 04 ) and 


A = 


/ 

0 

V 0 


0 =A0, 


o -i<3f-<s 

-zd> 0 
0 

+ ^ 0 


(3.4) 


0 \ 
0 


(3.5) 


It can be easily found that 

detA = (d > 2 + ^ 2 + ^ 2 ) 2 . 

The solution to the equation (13.41) can be written in the form 


0(0 = Texp Ai(T)drJ 0(fi), 


(3.6) 


(3.7) 


where 


A = 


/ 


V 


-iW 

0 

-i?¥ + y 

0 


0 

-i3> 

0 

i<3f + y 


-iw-y o 

0 iW-<3 

i® 0 


0 


i. 


(3.8) 


and 0 0i) is the solution at t — t \. As we have already shown, K — V 0 2 /V 2 for K = {/, I+ J, I — J} 
with trivial spinor-mass and K — Vq/V 2 for K = I for any spinor-mass. Since our Universe is 
expanding, the quantities S>, & and Sf become trivial at large t. Hence in case of K — I with non¬ 
trivial spinor-mass one can assume 0 (t \) = col (e _w Vi, e ~ tm ^h _ i _ e nn sv h ^ , w hereas for other 
cases with trivial spinor-mass we have 00 i) = col ( 00 , 0 ®, 0 ®, 0 °) with 0 ( ° being some constants. 
Here we have used the fact that <f> = m sp + The other way to solve the system (13.31) is given in 
[26], 


As far as equation for V, i.e., (12.471) is concerned, we solve it setting K = I as in this case we 
can use the mass term as well. 

Assuming 

F = = I>S 2 '^ (3.9) 

k k 

on account of S — Vq/V we find 


V = d>(y), q>(y)=xy 1/3 " 2/V3 + 


3 k 


m, 


sp 


Vb + 2£4(1 - n kK nk V 


2nk\/l—2iik 


2\ v (2K 3 -4/3) 


( 3 . 10 ) 


where X — 2 (m 2 — mn + n 2 ) X 0 

Let us now show the existence and uniqueness of the solution of the Eq. (13.101) . For this we 
study the right hand side of the Eq. (13.101), namely we check whether <£>0, V ) satisfies the Lipschitz 
condition. In doing so let us rewrite (13.101) as 


W = <J>(V), 
V — W. 


(3.11a) 

(3.11b) 
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Let (f*, V*) and (t*, IF*) be a particular pair of values assigned to the real variables (t, V ) and 
(f*, W*) such that within a rectangular domain D surrounding the point (z*. V*) and defined by the 
inequalities 

|f — f*| < A*, |V-V*|< j B*, (3.12) 

<J>(t, V) is a one valued continuous function of t and V. Indeed, for the other parameter fixed 
<f>(t, V) is a one valued continuous function. Recalling that V = 0 corresponds to a space-time sin¬ 
gularity and V is essentially non-negative, for any nontrivial value of V we conclude that <f>(/. V) \ 
has an upper bound M in D. We also define h = min (A*,B*/M). If h < A*, upon t we impose the 
additional condition \t — u\ < h. Let (t\, Vi) and (t 2 , V 2 ) are two points within D. 

The Lipschitz condition in this case 

s ][<&(*!, V\) - <J>(f 2 , V 2)] 2 + m - W 2 ) 2 <L^/(V 1 - V 2) 2 + (W! - W 2 ) 2 , (3.13) 

for any L > 1 follows from 

\! [<J>(fi, Vi) - <J>(f 2 , V 2)] 2 < Lsj (Vi - V 2 ) 2 . (3.14) 

Here L is a constant. Inserting <J>(t, V ) from (13.101) into the left hand side of (13.141) we find 


|4>(r 1 ,V 7 i)-4>(f 2 ,V 2 )| = 


x(v I/3-2N 3 _ v I/3-2N 3 


+ 3rc^X k (l-n k )V ( 


2 n k 


yl— 2 rik _yl— 2.11k 


(V 2 -V 1 ) 


li (y*)-( 2/3 +2i V 3) + £^(y 


*\- 2 n k 


(3.15) 


where V* G [Vi, V3], V) >0, Vi > 0. Here we also denote X\ — (1/3 — 2N^)X, and = 


2k(1 — nk ) (1 — 2rik)v } nk Here we used the mean value theorem. 


Since 


*i(y) 


r 0 

-(2/3+2N,) 


+LkK(y) 


- 2 n k 


val [Vi, V 2 ], such that 


is continuous, it possesses a maximum L in the inter¬ 


im, Vi) - <J>(f 2 , V 2 ) | < L\V 2 - V\ | (3.16) 

in the domain D. Further extending this study to other domains it can be shown that the condition 
13.16D holds in D = U(/ 1 >oy 2 <-[^i, V 2 ]. 

T hus w e see that <f>(t, V) is continuous and satisfies Lipschitz condition in the domain D. Hence 
Eq. (13.101) admits a unique continuous solution. 

Once the existence and u nique ness of the solution of (13.101) is proved, we can carry our study 
further. The first integral of (13.101) is 


V = <*>i(V), <J>i(V) = 


\ 


XiVW 3 -^) + 3 k 


m 


sp 


v 0 v+£ kkV^v^-”*)+c 


, (3.17) 


where we denote X\ — 6X/(4 — 6 ./V 3 ) and C is the constant of integration. The solution for V can 
be written in quadrature as 


dV 


X lV (4/3-2 n 3 ) +3k \ m v 0 V + £* X k V* lk V 2 ( 1 - n k)\ +C 


— t+tQ, 


(3.18) 
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with C and to being some arbitrary constants. 

In what follows we solve the Eqn. (13.101) numerically. In doing so we determine V (0) from 
(13.171) for the given value of V (0). 

To determine the character of the evolution, let us first study the asymptotic behavior of the 
equation (13.101) . It should be recalled that we have K — Vq/V 2 . Since all the physical quan¬ 
tities constructed from the spinor fields as well as the invariants of gravitational fields are in¬ 
verse function of V of some degree, it can be concluded that at any spacetime point where 
the volume scale becomes zero, it is a sin gular point [24]. So we assume at the beginning V 
was small but non-zero. Then from (13.101) we see that at t —> 0 the nonlinear term prevails if 
n k — n\\n\> max [1/2, A /3 + 1/3]. Recalling that we are considering an expanding Universe, at 
t —> 00 the volume scale should be quite large. In that case the nonlinear term prevails over the first 
term if n k — n 2 : n 2 < min [1/2, A /3 + 1/3]- For n k = no : 1 — 2 t 7 o = 0, i.e. no =1/2 the spinor field 
nonlinearity vanishes and the corresponding term becomes equivalent to the (effective) mass term. 

To define whether the model allows decelerated or accelerated mode of expansion we also plot 
deceleration parameter q defines as 

2 * < 3 - 19 ) 

which in view of (13.101) and (13.171) can be rewritten as 


q = 


y<F(y) 

*?(V) 


Xy4/3~ 2 N 3 + 3 k 

m sp V 0 V + 2I*4*(1 ~n k )vl nk V l 

—2n k 

XiVl 4 / 3 -^) +3 k 

m sp VoV + ^ k hV 2nk V 2 ^) 

+c 


(3.20) 


Now let us see what happens to deceleration parameter as t —* °°. As we have already estab¬ 
lished, for n 2 < 1/2 and 772 < 1/3 + A /3 the nonlinear tern prevails and in this case we find 


qm-(l-n 2 ) < 0 , 

whereas for N 3 < 1 /6 and nn > 1 /3 + A /3 we have 


(3.21) 


q 


X 


—( 2/3 — A/ 3 ) < 0 . 


(3.22) 


Thus we see that in both cases the Universe expands with acceleration. 

It should also be emphasized that for 771 > 1/2 and A /3 > 1/6 the mass term prevails asymptot¬ 
ically at t —y 00 and the Universe expands as a quadratic function of time, i.e., V | t 2 . 

The above analysis shows that the absence of mass term leads to constant deceleration param¬ 
eter, while for a time depending deceleration parameter the presence of a non-zero mass term is 

essential. _ 

Let us also see what happens to EoS (equation of state) parameter in this case. Inserting (13.91) 

into (12.23al) and (12.23bl) one finds the expressions for energy density £ = T 0 ° and pressure p = —A , 1 
(in this particular case as 7) 1 = A, 2 — T 2 )\ 

y y2nk y2 nk 

e=m sp y+£4^, p = £4(2n,-l)^. (3.23) 


In view of (13.231) for the EoS parameter we find 


y 2n k 

P _ 1 , 4 ( 277 ,- 1 )^ 
' >n s + 


(3.24) 


It can be shown that at the early stage of evolution the EoS parameter is dominated by the usual 
matter, while at later stage the dark energy becomes dominant. Moreover, at the absence of the 
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mass term the EoS parameter becomes a constant as in that case W — 2n k —\, whereas for a non¬ 
trivial mass term the EoS parameter is a variable function of time. Here we have exploited the fact 
that, at any concrete stage of evolution, one of the terms of the sum becomes predominant; hence 

others can be overlooked. _ 

There might be some question regarding the choice of nonlinearity in the form (13.91) . The 
reason lies on the fact that the spinor description of different kinds of fluid and dark energy such 
as ekpyrotic matter, dust, radiation, quintessence, Chaplygin gas, phantom matter etc. is in one 
form or the other is given by the power law of the invariants of spinor field. While t he sp inor 
description of fluid or dark energy leads to the elimination of mass term, the choice (13.91) still 
allows us to study the role of spinor mass on the evolution of the Universe. To show this let us 
recall that only in case of K = I — S 2 we could express K in terms of V with a non-trivial mass 
term in the Lagrangian. So setting F = F(S) from 


W = 


in view of (I2.23al) and (I2.23bl) one finds [43-461 


V f 


i+w 


Vo 


F — AS 1 " 1 "™ 1 — mS — X °. ... — m 

yl+W y 


(3.25) 


(3.26) 


that corresponds to dust (W = 0), radiation (W = 1/3)), hard Universe (W G (1/3, 1)), stiff 
matter (W — 1), quintessence (VE G ( — 1/3, —1)), cosmological constant (VE = — 1 ), phan¬ 
tom matter (VE < —1), and ekpyrotic matter (XV > 1), respectively. Inserting (13.261) into (12.71) 
one finds that the mass term in this case vanishes, while the spinor field nonlinearity given by 
(13.91) does not. In this case for energy density and pressure we find £ = XVq +w /V l+W and 
p = AWV 0 1+m/ /V' 1+v ' / , respectively. EoS parameter in this case is a constant by definition, while 
in absen ce o f the mass term the deceleration parameter also becomes a constan t. N evertheless one 
can use (13.91) with the trivial mass term in the Lagrangian and the sum (£/,) in (13.91) can be viewed 
as multi-component source field with k standing for different types of matter and dark energy such 
as ekpyrotic ma tter, dust, r adiatio n, quintessence, Chaplygin gas, phantom matter etc. 

Comparing (13.91) with (13.261) one finds 2 n k — W + 1. Further setting the value of XV for 
different fluid and dark energy we find the corresponding value of ny. dust (n k = 1/2), ra¬ 
diation (n k = 2/3)), hard Universe (n k G (2/3, 1)), stiff matter (n k = 1), quintessence 
(n k G (0, 1/3)), cosmological constant (n k — 0), phantom matter (n k < 0), and ekpyrotic mat¬ 
ter (n k > 1), respectively. It was shown earlier, when n k — 1/2 the corresponding term can be 
added to the mass term. So we can conclude that the term with n k — 1/2 which also describes dust 
behaves like a mass term. 

One of the principal advantage of using spinor description of source field lies on the fact that in 
this case one needs not think about whether two or more components considered can be separated. 
To show that let us write the Biacnhi identity that leads to 

y.v = ry+r; v y - = o. (3,27) 

which for the metric (12.41) on account of the components of the energy-momentum tensor takes the 
from 

£ + — (£ + p) — 0. (3.28) 

Inserting £ and p from (I2.23al) and (12.23bl) from (13.281) one finds 

m sn d , , Fv d , i. 

In case of K = I — S 2 (13.291) fulfills identically thanks to (12.491) . i.e., SV — const, and KV 2 = 
const., whereas in the case when K takes one of the following expressions {J, I + J, I — J}, for a 
massless spinor field (13.291) fulfills identically thanks to (12.511) . (12.541) and (12.581) . i.e., KV 2 — const. 
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Hence if we use spinor description of different fluid and dark energy simulated from corresponding 
equation of sate, the Bianchi identity will be fulfilled identically without invoking any additional 
condition. 

To this end let us solve the equation for V numerically. For simplicity let us consider system 
with two components only. In this case we have 


V = <*>(V), 


(3.30) 


3 /c* 

<F(V) = xv l/3 ~ 2N3 + — \m sp V 0 + 2?n(l-n l )VQ ni V 1 - 2ni + 2A 2 (1 -n 2 )V 0 2 " 2 V 1_2M2 

with the first integral 

V = <J>t(V), (3.31) 

d>j(y) = i i y( 4 /3- 2iV 3)+3K : L S pVoV+AiV 0 2wi y 2 ( 1 - ni )+>i 2 v r 0 2M2 y 2 (i“ M 2) + cl, 


Since we are interested in qualitative picture here, so we set the value of problem parameters 
very simple. In doing so we set the value of N$, n\ and « 2 in such as way that none of the four terms 
in the right hand side of (13.301) merge with the others. Beside this we will consider the coupling 
constants X\ and A 2 with different signs. The i nitial value V(0) is taken to be small but non-zero 
in such a way that the right hand side of (13.311) remains non-negative. For the given initial value 
V (0) is defined from (13.311) . 

From (13.311) it can be easily established that only in case when both Ai and A? are positive the 
model allows ever expanding solution, whereas, if one of the A/’s is negative, the non-negativity of 
the expressions under the square-root imposes some restrictions on the value of V. A negative Ai 
generates the minimums while the negative A 2 generates the maximums. In case if the minimum 
occurs at a negative value of V we have a Universe that expands to some maximum value and then 
contracts before ending in a Big Crunch. If both maximums and minimums are non-zero, we have 
periodic solution with no beginning and no end. In both cases A 2 <0. If A 2 > 0 independent to 
whether Aj is positive or negative, we have ever expanding solution. 

For simplicity we set m = 1, n = 2, Xo = 1, Vo = 1, w sp = 1, Co = 1, fc = 1. Fixing N3 = —1/3 
from n\ > max [1/2, N 3 + 1/3] we set 11 \ =3/2 (ekpyrotic matter) and n\ = 2/3 (radiation), while 
from n 2 < min [1/2, N 3 + 1/3] we set n 2 = 1/4 (quintessence) and n 2 = —1 (phantom matter). 
As far as coupling constants are concerned we consider two cases with A/ = {1,-0.001} and 
A 2 = {1, — 1}. The initial value of V (0) is taken to be V (0) = 0.01. 

In Fig. Q] we have illustrated the evolution of volume scale V for the Universe filled with 
massive spinor field with nj = 2/3, n 2 = 1/4, Aj = 1 and A 2 = 1. We draw the picture of evolution 
of V in the Figs. [2] [3] and [4] for {n\ = 2/3, « 2 = — 1, A/ = 1, A 2 = —1}; {n\ — 3/2, « 2 = 1/4, A/ = 
—0.001, A 2 = 1}, and {n\ = 2/3, n 2 = — 1, Ai = —0.001, A 2 = —1}, respectively. In Figs. [5] and [6] 
we have illustrated the evolution of deceleration parameter for positive A 2 only. 


IV. CONCLUSION 


Within the scope of Bianchi type-VI spacetime we study the role of spinor field on the evolution 
of the Universe. It is found that in this case the non-diagonal components of the energy-momentum 
tensor of spinor field, unlike in the cases Bianchi type I [48] and Bianchi type-V/o [491, does not 
lead to the elimination of spinor field nonlinearity and the mass term in spinor field Lagrangian. 
Depending of the sign of self-coupling constant the model in this case allows eilher late time 
acceleration or oscillatory mode of evolution. 
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FIG. 1: Evolution of the Universe filled with massive spinor field with n\ = 2/3, n 2 = 1/4, Ai = 1 and 
A 2 = 1 
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FIG. 2: Evolution of the Universe filled with massive spinor field with n\ = 2/3, nj = —1, A| = I and 
A2 = — 1 . 
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FIG. 4: Evolution of the Universe filled with massive spinor field with n\ = 2/3, n 2 = —1, Ai = —0.001 
and A 2 = — 1. 
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FIG. 5: Plot of deceleration parameter q with n\ = 2/3, m = 1/4, Ai = 1 and A 2 = 1 
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FIG. 6: Plot of deceleration parameter q with n\ = 3/2, «2 = 1 /4, Aj = —0.001 and 2,2 = 1. 





